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Abstract. There are two types of 3D shears in channel flows: (U(y, z), 0, 0) 
and (U(y), 0, W(y)). Both are important in organizing the phase space struc- 
tures of the channel flows. Stability criteria of the type I 3D shears were 
studied in [Li, 2010]. Here we study the stability criteria of the type II 3D 
shears. We also provide more support to the idea of resolution of a turbulence 
paradox, introduced in [Li and Lin, 2010], by studying a sequence of type II 
3D shears. 



1. Introduction 

Studying the phase space structures of channel flows is an important and emerg- 
ing area. Most of the existing works in this area is numerical. Like every other 
dynamical system study in a phase space, non-wandering objects like fixed points, 
periodic orbits etc. play a fundamental role in organizing the phase space structure. 
When the Reynolds number is infinite (i.e. zero viscosity), the corresponding 3D 
Euler equations have two types of steady shears as fixed points in the phase space. 

• Type I 3D shears (U(y, z), 0, 0), 

• Type II 3D shears (U(y),0, W{y)); 

where the boundaries of the channel are in the y-direction. When the Reynolds 
number is not infinite but large, these shears turn into slowing drifting states under 
the corresponding 3D Navier-Stokes dynamics. By the concepts of rate condition 
and normal hyperbolicity of Fenichel pQ [2], these slowing drifting states can be 
crucial in the transition to turbulence. Also the 3D Navier-Stokes dynamics has 
fixed points like the so-called lower and upper branches [9]. When the Reynolds 
number approaches infinity, the lower branch approaches one of the type I 3D 
shears. How to distinguish this particular one from the rest of the type I 3D shears 
is an interesting question (also posted in a list of problems by Yudovich |10j ). A 
condition satisfied by this particular type I 3D shear was derived in [7]. Again 
like every other dynamical system study in a phase space, the stability of these 3D 
shears is crucial in understanding the phase space structure. The stability criteria 
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for type I 3D shears were studied in [5]. Here we shall study the stability criteria 
for type II 3D shears. It turns out that the linearized 3D Euler equations can be 
casted into formally the same form as the classical Rayleigh equation for 2D shears 
(U(y),Q). But the nature of the stability criteria is fundamentally different from 
that for the 2D shears. One can ask the question: What is the "percentage" among 
e.g. all type I 3D shears, that is unstable? The author's conjecture is: 

• The unstable percentage of type II 3D shears > the unstable percentage 
of type I 3D shears > the unstable percentage of 2D shears. 
Finally, we shall provide more support to the idea of resolution of the turbulence 
paradox, introduced in [6] by studying a sequence of type II 3D shears. The turbu- 
lence paradox is also called Sommerfeld paradox which roughly says that the linear 
shear in the plane Couette flow is linearly stable for all values of the Reynolds 
number, while in experiments, transition from the linear shear to turbulence occurs 
when the Reynolds number is large enough. For more details on the turbulence 
paradox and our resolution, see Section |4] and the paper [6] . 

2. Necessary Conditions For Instability 

The inviscid channel flow is governed by the 3D Euler equations 

(2.1) d t Ui + UjUij = -p t i, u iti = 0; 

where {u\,U2,u 3 ) are the three components of the fluid velocity along (x,y,z) 
directions, and p is the pressure. The boundary condition is the so-called non- 
penetrating condition 

(2.2) U2(x, a, z) = 0, U2(x, b, z) = 0; 

where a < b are the boundary locations of the channel in y-direction. 

We start with the type II 3D steady shear solutions of the 3D Euler equations: 

ui = U{y), u 2 = 0, u 3 — W(y), p — po (a constant). 

Of particular importance are those profiles which also satisfy the non-slip boundary 
condition 

U(a) = a, U(b) = /3, W(a) = W(b) = 0; 
where a and (3 are the velocities of the two walls of the channel. Such profiles may 
be the viscous limiting profiles when the viscosity approaches zero. Linearize the 
3D Euler equations with the notations 

Ul = U(y) + [e iklx+ik3Z - i(7t u(y) + c.c] , u 2 = e lklX+lk3Z - l<Tt v(y) + c.c, 

u 3 = W(y) + [e iklx+ik3Z - i(7t w{y) + c.c] , p -> Po + [^ x+ik »*- iat p(y) + c.c] ; 

where h\ and k 3 are the wave numbers, and a is a complex constant; we obtain the 
linearized 3D Euler equations 



(2.3) i(kiU + k 3 W-a)u + U'v = -ik 1 p, 

(2.4) i(kiU + k 3 W - a)v = -p' , 

(2.5) i{k t U + k 3 W - a)w + W'v = -ik 3 p, 

(2.6) ik\u + v + ik 3 w = 0. 



Two forms of simplified systems can be derived: 

fci U" + kuW" 
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with the boundary condition v(a) — v(b) = 0; and 

(2.8) {fall + k 3 W - a) 2 [(fall + k 3 W - a)- 2 p']' = {kf + k 2 )p, 

with the boundary condition p'(a) — p'(b) = 0. 

Notice that formally equation (12 . T[) is in the same form as the clasical Rayleigh 
equation where the 2D shear U(y) is now replaced by k\U + k^W . We have the 
following extension of the classical Rayleigh's inflection-point theorem. 

Theorem 2.1. // (ki,ks) is an unstable mode, then 

(2.9) (kiU" + k 3 W") (y») = for some y« e (a, b). 



Remark 2.2. Even though the Rayleigh's inflection-point theorem is a special case 
of the above theorem, the claim in the above theorem in general is fundamentally 
different from that of Rayleigh. Rayleigh's claim basically says that if (U(y),0) 
is linearly unstable, then U{y) has an inflection point. Here on the other hand, 
condition (|2.9[) can be satisfied by most of type II 3D shears {U(y), 0, W(y)). In 
fact, one can choose k\ — W"(y*), k 3 = — U"(y*) for any y* 6 (a, b), then (|2.9[) 
is satisfied. Therefore, in general (I2.9|) is a very weak necessary condition for the 
linear instability of (U(y),0, W(y)). This may indicate that type II 3D shears are 
more often to be unstable than 2D shears. The author's conjecture is that type I 
3D shears are inbetween in terms of frequency of instability. 

Proof. Multiplying (|2~7l) by v, 

("0) / IWf + « + *»h>fl * + [ XffiV * " °' 

the imaginary part of which is 

( 2 - U ) ^ / 17 rr I 7 w 12 M d y = ' 



□ 



/cif/ + fc 3 VF-cr| 2 

where a = a r + icj . Thus 

(fci(7" + k 3 W") (y*) = for some € (a, 6). 



We also have the following extension of the Fj</>rtoft's theorem. 
Theorem 2.3. If (k\,k 3 ) is an unstable mode, then 

G"(yo) [G(y ) ~ G(y*)] < for some y £ (a, b), 
where G — k\U + k 3 W and y* is the point at which G" = given by Theorem \2.1 

PROOF. The real part of equation (|2.10p is 

f a G "\G-J^ V?dv = " jf [ |wf + {kl + ^ 2)H2 ] dy - 
From (|2.1ip . one has 

G 



[oy - G(y*)] f 



lG _ al M dy = o, 
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where is chosen to be the point given by Theorem 12.11 G"(y*) = 0. The above 
two equations imply that 

b G"(y)[G(y)-G^)} t 2 

\G{y)-o\* H dy< °' 

and the theorem is proved. □ 

Next we show the extension of Howard's semi-circle theorem. 

Theorem 2.4. If (ki,k 3 ) is an unstable mode, then its corresponding unstable 
eigenvalue lies inside the semi- circle in the complex plane: 

M + m\ 2 , [M-m\ 2 

* — —) ' 

where again a = a r + iui, M = maxj, e [ a y (k\U + k 3 W), and m — min^gu b -\{k\U + 
k 3 W). 

Proof. Multiply (|2.8[) with p(G — er)~ 2 , integrate by parts, and split into real 
and imaginary parts; we obtain that 

r-b r-b 

(2.12) / GQdy = a r Qdy, 

J a J a 

pb pb pb 

I G 2 Qdy = 2a r GQdy + (of - a 2 ) Qdy 

J a J a J a 

f b 

(2-13) ={o 2 r +° 2 ) \ Qdy 

J a 

by (|2~T2|) . where 

G = kiU + k 3 W, 

Q = \G-o-\-*[\pf + (kl + kl)\pf]. 



Let 



then 



M = max G, m — min G, 

y£[a,b] y&[a,b] 



b 

(G - m)(M - G)Qdy > 0. 



Expand this inequality and utilize (|2.12[) - (|2.13[) . we arrive at the semi-circle in- 
equality in the theorem. □ 

3. Sufficient Conditions For Instability 

The zeroth mode (ki, &3)=(0, 0) is trivially neutrally stable, so our interest 
is focused upon non-zero modes. Without loss of generality, we assume k\ ^= 0. 
Equation (|2.7p can be re-written in the form 



U" + ^-W" 

V" r-^ V = k\ 



1 



We will keep k = ks/ki fixed and vary k\. Denote by 

H = U + kW. c = a = kiy/i + k 2 , 
ki 
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we obtain the following equivalent form of (|2.7|) . 

H" 

(3.1) v" - — v = a 2 v, 

H — c 

with the boundary condition v(a) — v(b) = 0, which is in the same form as the clas- 
sical Rayleigh equation 6 . Thus we have the extension of the Tollmien's theorem 
on a sufficient condition for instability. 

Theorem 3.1. If H'(y) ^ for all y € (a, b), H"(y„) = for some y„ € (a, b), 
and the Sturm- Liouville operator 

Lv = —v 



H - H{y* 



has a negative eigenvalue under the Dirichlet boundary condition v(a) — v(b) = 0; 
then the equation \3.1\) has an unstable eigenvalue (in fact, an unstable curve c = 
c(a) for a in some interval). 



4. Turbulence Paradox 

Turbulence paradox is also called Sommcrfcld paradox. It originated from 
Sommerfeld's analysis which concluded that the linear shear in plane Couette flow is 
linearly stable for all values of the Reynolds number; whereas in fluid experiments, 
perturbations of the linear shear often lead to transition to turbulence. Such a 
paradox is universal among fluid flows, e.g. pipe Poiseuille flow, plane Poiseuille 
flow etc. . A resolution of this paradox is given in [6] . The main idea of the resolution 
is to show that even though the linear shear is linearly stable, states arbitrarily close 
to the linear shear can still be linearly unstable. Here different norms are crucial. 
It is shown in [6] that the sequence of 2D shears (U n (y),0) is linearly unstable for 
all n and large enough Reynolds number (including infinity), where 

(4.1) U n (y) =y+^ sin(4mry), Q-L < A < i-) . 

Here U n (y) approaches the linear shear U = y in the L°° norm of velocity, but not 
in the L°° norm of vorticity. Notice also that when the Reynolds number is not 
infinite but large, the shears (U n (y),0) are not steady states rather slowly drifting 
states. The linear instability mentioned above is predicted by the Orr-Sommerfeld 
operator (linearized 2D Navier-Stokes operator) when the shears (U n (y),0) are 
viewed frozen. Such an instability is in the spirit of Fenichel's rate condition and 
normal hyperbolicity pQ [2] . Such an instability is also observed numerically [4] . 

In this paper, we would like to add more support to the idea of resolution 
mentioned above by considering the type II 3D shears (U(y),0, W{y)). 

The viscous channel flow is governed by the Navier-Stokes equations 

(4.2) d t Ui + UjUij = —p ti + euijj, u»,j = 0; 

where again [u\, U2, U3) are the three components of the fluid velocity along (x, y, z) 
directions, p is the pressure, and e = 1/R is the inverse of the Reynolds number R. 
The boundary condition is 



(4.3) u\ (x, a, z) = a, u±(x, b, z) = /3, Uj(x, a, z) = Uj(x, b, z) = 0, (J = 2, 3); 
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where a < b, a < (3. For the viscous channel flow, the type II 3D shears mentioned 
above are no longer fixed points, instead they drift slowly in time (sometimes called 
quasi-steady solutions): 



[e ad lu{y),Q^ td lW{y)) 



By ignoring the slow drift and pretending they are still fixed points (or by using 
artificial body forces to stop the drifting), their unstable eigenvalues will lead to 
transient nonlinear growths as shown numerically [4]. A better explanation here 
is to use the theory of geometric singular perturbation pQ [2] [3] (§]. The slowly 
drifting 3D shears altogether form a locally invariant slow (center) manifold. The 
normal direction growth rate (or decay rate) of this slow manifold has a persistence 
property (i.e. robust). Thus the growth rate can be estimated by ignoring the slow 
drift. The geometric singular perturbation theory implies the transient nonlinear 
growth induced by the linear growth rate. 

The corresponding linear Navier-Stokes operator at (U {y),0, W(y)) is given by 
the following counterpart of (|2.3p - (|2.6|) . 

(4.4) i(kiU + k 3 W - a)u + U'v = -ik lP + e[u" - {k\ + fcf)u], 

(4.5) i(k x V + k 3 W - a)v = -p' + e[v" - (fcf + k\)v}, 



3) 

(4.6) i(kiU + k 3 W - a)w + W'v = -ik 3 p + e[w" - (k 2 + fcf )iu] 

(4.7) ik x u + v' + ik 3 w = 0. 

The simplified system as the counterpart of (|2.7p is 



(4.8) - 

la 



2 

dy 2 



2 r d 2 

v + H"v-(H-c) 



dy 



2~ a 



v = 0, 



with the boundary condition v — v' — at y — a, b; where e = eyl + n 2 , as 
before k = k 3 fk\ (fixed), a = k\ yl + k 2 , c = cr/k\, and H = U + kW. Now both 
equation (|3.1j) and equation (|4. 8[) are formally in the same form as those in [6] . We 
can specify 

a = 0, 6=1, a = 0, P = l\ 
and consider the following sequence of type II 3D shears 

(4.9) (U n (y),0, W n (y)) = (y, 0, - sin(4wry) 

The problem of linear instability of the type II 3D shears (|4.9p is then casted into 
the same problem as that of the 2D shears (|4.1j) . By the results of [6], we have 

Theorem 4.1. For any A > and any integer n > 1, the type II 3D shears 
\4-9$ are linearly unstable under the 3D Euler dynamics. Specifically, there exists 
a 2D unstable eigenmode surface (ki,k 3 ;a(ki,k 3 )) with Im{a{k\, k 3 )} > for the 
equation \2. 7| ), stemming from neutral modes of the form (ki(n),k 3 (n); ^ki(n)) 
where [fci(n)] 2 + [k 3 (n)] 2 > Cn 2 (C > is independent of n). The corresponding 
eigenf unctions are in C°°(0, 1). 



Proof. Fix k = ks/ki such that 

1 1 
24^ 

then the problem is reduced to that of Theorem 3.2 in [6]. □ 



11 1 

< kA < — , 

2 4?r 4tt 
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Theorem 4.2. Let (k®,k®;o ) be a point on the 2D unstable eigenmode surface 
given by Theorem \4-l\ Then when e is sufficiently small, there exists an unstable 
eigenmode (ki,k®;cr*) with Im{<7*} > for the equation with H given by 

(fpp. When e -> 0+ , a* -> a . 



Proof. For the fixed (fej, k®), K has a fixed value; then e and e are equivalent. 
The problem is reduced to that of Theorem 4.1 in [6]. □ 

Remark 4.3. Notice that the sequence of type II 3D shears (|4.9p are linearly 
unstable for all A > 0, while the sequence of 2D shears (|4.ip are proved linearly 
unstable for i -X- < A < The sequence of type II 3D shears (|4.9|) also approaches 
the linear shear U = y in the L°° norm of velocity, but not in the L°° norm of 
vorticity. 
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